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Abstract. We present a grid of evolutionary tracks for low- 
mass white dwarfs with helium cores in the mass range from 
0.179 to 0.414 Mq. The lower mass limit is well-suited for 
comparison with white dwarf companions of millisecond pul- 
sars. The tracks are based on a 1 Mq model sequence extending 
from the pre-main sequence stage up to the tip of the red-giant 
branch. Applying large mass loss rates at appropriate positions 
forced the models to move off the giant branch. The further 
evolution was then followed across the Hertzsprung-Russell di- 
agram and down the cooling branch. At maximum effective 
temperature the envelope masses above the helium cores in- 
crease from 0.6 to 5.4 • 10 -3 M for decreasing mass. We 
carefully checked for the occurrence of thermal instabilities of 
the hydrogen shell by adjusting the computational time steps 
accordingly. Hydrogen flashes have been found to take place 
only in the mass interval 0.21 <> M/M© <, 0.3. 

The models show that hydrogen shell burning contributes 
significantly to the luminosity budget of white dwarfs with he- 
lium cores. For very low masses the hydrogen shell luminosity 
remains to be dominant even down to effective temperatures 
well below 10 000 K. Accordingly, the corresponding cooling 
ages are significantly larger than those gained from model cal- 
culations which neglect nuclear burning or the white dwarf pro- 
genitor evolution. 

Using the atmospheric parameters of the white dwarf in the 
PSR J1012+5307 system we determined a mass of M = 0.19± 
0.02 Mq and a cooling age of 6 ± lGyr, in good agreement 
with the spin-down age, 7 Gyr, of the pulsar. 

Key words: Stars: evolution - Stars: interiors - White dwarfs 
- Binaries: general - Pulsars: PSR J1012+5307 



1. Introduction 

The theory of stellar evolution shows that white dwarfs (WD) 
with a carbon-oxygen core (CO-WDs) can be well understood 
in terms of the evolution of single stars with initial masses be- 
tween 1 and about 6 ... 8 M©. After central hydrogen burning, 
stars in this mass range develop first a helium core surrounded 
by a hydrogen-burning shell (red-giant branch or RGB) and 
later, after completion of central helium burning, a compact 
carbon-oxygen core surrounded by two shells burning hy- 
drogen and helium, respectively (asymptotic giant branch or 
AGB). At high luminosities strong stellar winds erode the stel- 
lar envelope effectively. When the envelope mass drops below 
a few 10~ 2 M© the model leaves the AGB, becomes a central 
star of a planetary nebula and finally cools down as a white 
dwarf. The final masses range typically between about 0.5 and 
1 M (see e.g. Iben &Renzini|1983fc Sch onberner |1 979] |1983 
Vassiliadis & Wood |l993L 119941; Blocker [l995aL |1995bh. 



White dwarfs with smaller masses cannot be produced by 
single-star evolution since their progenitors would have ini- 
tial masses below 0.5 M Q and thus do not finish their main 
sequence evolution within a Hubble time. Instead one has to 
invoke mass transfer in a close binary system whereby the 
donor's evolution towards central helium ignition is choked off 
the so-called case B mass exchange (Kippenhahn & Weigert 



1967). The remnant is a low-mass object (M < 0.5 M Q ), con- 



sisting of a helium core with a hydrogen-rich envelope still 
burning hydrogen at its bottom, which contracts slowly towards 
a white dwarf configuration, quite similar to the more massive 
and more luminous remnants from the AGB. 

Two known common types of binary systems can con- 
tain such helium white dwarfs (He-WDs). Firstly the so-called 
double-degenerate systems, where both components are white 
dwarfs (see e.g. Saffer et al. 1998). Examples for such systems 
with one or even two He-WDs are WD 0957-666 (Bragaglia et 
al. Il990|; Moran et al. Il997|), WD 01 35-052 (= L870-2 = EG1 1 , 



Send offprint requests to: T. Driebe 

* Tables 5-11 are only available in electronic form at the CDS via 
ftp 130.79.128.5 



Saffer et al. 1988), PG 1101+364 (= Ton 1323, Marsh 1995) 



and the five systems studied by Marsh et al. ( |1995| ). The sec- 
ond type are the so-called millisecond pulsar systems (hereafter 
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MSP), where the He-WD is the companion of a pulsar with a 
high rotational speed in a nearly circular orbit. 

So far about three dozens of these systems are known. Most 
of them have a low-mass white dwarf companion (see for ex- 
ample Bail es & L orimer 1995, Camilo et al. 1996, Ray et al. 
1996, Lyne 1996 ). Reviews of MSP properties are given, e.g., 
by Ver bunt ( |1993| ), Shor e et a l. (1994), Phinn ey & Kulkarni 
d l994 |), Rappap ort et a l. ( |1995| ), Nicastro et al. ( |1995| ), Camilo 
( 1996 ) or Lyne ( 1995 ). The evolution of the progenitor system 
can be explained by highly non-conservative mass exchange 
events (common envelope evolution). A detailed description of 
the binary scenario which leads to the formation of a MSP sys- 
tem has recently been given, e.g., by Tauris & Bailes (1996), 
Tauris (1996) or Ergma & Sarna (1996). In progenitor systems 
to MSPs the more massive component (Mi « 10 M Q ) fills its 
Roche volume, looses a substantial amount of mass but follows 
the evolution towards a supernova explosion with the formation 
of a neutron star which might be detected as a radio pulsar. The 
less massive secondary component (M2 ~ 1 M©) continues its 
evolution, evolves off the main sequence to become a red giant 
and will eventually also exceed its Roche lobe to spill its mass 
onto the neutron star. This accretion of matter spins up the neu- 
tron star to rotational periods of the order of milliseconds. In 
this phase the system can be identified as a low mass X-ray bi- 
na ry syst em (LM XB) lik e Cyg X-2 or 2S 0921-63 (Webbink et 



al. 1983; Verbunt 1993; for a population synthesis of LMXBs 



see Iben et al. 1995 ). Continuing mass exchange of the sec- 
ondary causes its envelope mass to fall below a critical value. 
The system becomes detached again, and the secondary evolves 
to a white dwarf with a helium core while the pulsar's rotation 
slows down again by energy loss at the expense of rotational 
energy and the pulsar properties are activated again. 

Using smaller initial mass ratios the scenarios for the emer- 
gence of MSs also describe, in principle, the formation of 
double degenerate systems, i.e. systems containing two white 
dwarfs (CO+CO, He+CO, He+He). For a discussion of these 
systems see f or example Webbink ( 1984 ), Iben & Tutukov 
( |l984a[ , |1986|), C ameron & Iben ( |1986| ) and, more recently, 
Sarna et al. (h 996|). 



MSP systems open the unique possibility to check the spin- 
down theory of pulsars by independent age determinations of 
the white dwarf components once the atmospheric parameters 
of the latter are known. PSR J1012+5307 is so far the only sys- 
tem where the white dwarf has been studied with sufficient pre- 
cision as to allow for comp aring spin-down and cooling ages 
(see e.g. Lorimer et al. 1995 ) and to determine the components' 
masses. 



PSR J1012+5307 was identified by Nicastro et al. ( |1995[ ). 
They found a rotational period P rot = 5.26 ms and a com- 
pact companion with a minimum mass of 0. 1 1 M© in an orbit 
with period P or b = 14.5 h. The surface parameters of the com- 
panion have been determined by van Kerkwijk et al. (1990 
and Callanan et al. (1998). Due to different atmospheric mod- 
els the derived parameters differ somewhat: van Kerkwijk et al. 
( |1996| ) found log.g = 6.75 ± 0.07 and T cff = 8550 ± 25 K 
while Callanan et al. <\l99$) give log g = 6.34 ± 0.20 and 



T e ff = 8670 ± 300 K. According to these atmospheric param- 
eters the companion is a low-mass helium white dwarf as sug- 
gested by the theoretical scenarios sketched above. Van Kerk- 
wijk et al. ( |1996| ) estimated a mass of 0.16 ± 0.02 M Q by 
extrapolating from more massive carbon-oxygen models. The 
same result was found by Callanan et al. (1998 ) from their anal- 
ysis and white dwarf modelling. 



Lorimer et al. (1995) determined the spin-down age r = 
7.0 Gyr assuming Pq <C P I0 t, where Pq is the initial pulsar 
period after the end of the spin-up phase. The age of the MSP 
can be estimated as well by the cooling age of the WD com- 
panion. However, the estimates based on existing white dwarf 
calculations are controversial: while Lorimer et al. (1995) de- 



rived a cooling age for the white dwarf of about 0.3 Gyr, which 
is a factor 20 lower than the pulsar's spin-down age, Alberts 



et al. (1990 give a value of 7 Gyr, in full agreement with the 



spin-down age of the pulsar. 

To solve the problem of age determination of MSPs two 
different strategies are possible: The first concerns the critical 
examination of the spin-up and spin-down theories for pulsars. 
On one hand, an uncertainty in the spin-down age determina- 
tion arises from the assumption Pq <C P- Violations of this 
assumption result in lower spin-down ages (see for example 
Camilo et al. 1994). On the other hand, Burderi et al. ( |1996j ) 
showed that by assuming an accretion-induced field decay in- 
stead of spontaneous field decay substantially lower spin-down 
ages can be obtained. Following their conclusions MSP sys- 
tems would thus have characteristic ages of a few 10 8 years 
only. In the particular case of PSR J1012+5307 Burderi et al. 



(1996 ) found a spin-down age of about 3.7 • 10 yr. 



The other strategy for deriving corresponding ages involves 
the reexamination of the evolutionary models of low-mass 
white dwarfs. Studies of such objects are rare, in particular for 
M < 0.2 M Q as is appropriate fo r PSR J1012+5307. Sev- 
eral authors, i.e. Kippenhahn et al. (1967), Kippenhahn et al. 
d!968| ), Ref sdal & Weigert ( |1969| ), Gian onne e t al. ( |1970| ), Iben 
& Tutukov d!986[ ) and Castellani et al. ( |1994| ) computed mod- 
els of He-WDs with M > 0.2 M by simulating mass ex- 
change in close binaries. Typical is that all of these studies have 
encountered thermal instabilities (or hydrogen flashes) in mod- 
els with M > 0.25 M Q due to unstable hydrogen burning on 
the cooling branch. 



Chin & Stothers (1971) and Webbink (1975) constructed 



He-WD models by following the evolution of single stars of 
the appropriate low masses. Webbink (1975) calculated a grid 
of sequences in the mass range 0.1 <C A^f/M© <C 0.5. However, 
due to large time steps used his tracks do not show any thermal 
instabilities. From Webbinks calculations one can infer that the 
white dwarf companion in the PSR J1012+5307 system should 
have an age of r w 10 Gyr. Chin & Stothers (1971) did not 
include hydrogen shell burning in their white dwarf models. 
According to their M = 0.2 M© track, the white dwarf in the 
PSR J1012+5307 system is only 0.3 Gyr old. Hereafter, models 
with the explicit assumption L nuc = will be referred to as 
contraction models. 
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Alberts et al. (1996) calculated He-WD models between 
0.17 < M/M© < 0.25 by simulating binary evolution with 
low-mass components. They found no thermal instabilities of 
the hydrogen burning shell. Their models give for the PSR 
J1012+5307 companion's mass and age M = 0.185 M Q and 
t = 7 Gyr, respectively. 



Althaus & Benvenuto (1997) and Benvenuto & Althaus 



( 1998 ) presented models in the mass range 0. 15 < Mj M < 
0.5. Althaus & Benvenuto ( |1997[ ) treated He-WDs without any 
hydrogen envelopes, and thus did not allow for nuclear burning. 
Benvenuto & Althaus ( |1998[ ) and Hansen & Phinney ( |1998^ ) 
(0.1 < M/Mq < 0.5) started their calculations considering 
hydrogen burning, but found it to be insignificant. All sets of 
tracks of these authors display very similar cooling properties, 
and one can estimate a cooling age of about r < 0.5 Gyr for the 
white dwarf in the PSR J1012+5307 system. Note, that, how- 
ever, their initial models differ from those based on calculations 
of RGB progenitors with mass loss. 

It is obvious from this comparison that the cooling prop- 
erties of low-mass white dwarfs are extremely sensitive to the 
initial conditions. We suggest that initial models closest to real 
He-WD progenitors result from explicit modeling of their evo- 
lutionary history which determines the model's envelope mass 
and thermomechanical structure. Therefore we have calculated 
an extensive grid of model sequences with full consideration 
of nuclear burning. Our investigations revealed that the simul- 
taneous consideration of these two aspects, i.e. nuclear burn- 
ing and the evolutionary history, results in considerably longer 
evolutionary ages of the white dwarfs. Consequences for the 
mass determination are evident as well. A full presentation of 
this evolutionary grid, with an extensive discussion of thermal 
instabilities and the mass-radius relation is deferred to a forth- 
coming paper. Here we will discuss only those model proper- 
ties that are important to interpret millisecond pulsar systems. 

This paper is organized as follows: in Sect. || we describe 
the stellar evolution code and the method used to calculate evo- 
lutionary models presented in this paper. Results are given in 
Sect. and applied to the particular system PSR J1012+5307. 
Sect. |[ summarizes our results. 

2. The evolutionary calculations 

The evol utionar y code used is essentially the one described by 
Blocker (1995a), with several modifications. Nuclear burning 
is accounted for via a nucleosynthesis network including 30 
isotopes with all important reactions up to carbon burning sim- 
ilar as in El Eid (1994). The most recent radiative opacities by 
Iglesias et al. ( |1992[ ) and Iglesias & Rogers ( |1996[) [O PAL], 
supplemented by those of Alexander & Ferguson ( 1994) in the 
low-temperature region, are employed. Diffusion is not con- 
sidered. The initial composition is (Y, Z) = (0.28, 0.02), the 
mixing length parameter a = 1.7 follows from calibrating a 
solar model. The Coulomb corrections to the equation of state 
are those given by Slattery et al. (1982). We note in passing 
that for the results presented here the quality of the Coulomb 
correction treatment is of no significance. 



For comparison we have also computed evolutionary se- 
quences based on different input physics: 

- sequences using the older opacities of Cox & Stewart 
dl965ai|l965bl|l970| ) [CS]; and 

- sequences using equilibrium nuclear reaction rates for hy- 
drogen and helium burning. 

In the following we will only refer to the sequences based on 
the OPAL opacities with full employment of the nuclear net- 
work if not mentioned differently. 

2.1. Method of calculation 

Because we primarily focused our study on the cooling be- 
haviour of low-mass white dwarfs and the implications for the 
mass-radius relation we did not calculate the mass exchange 
phases duri ng the RGB evolution i n deta il (see for this subject, 
e.g., Tauris 1996, Tauris & Bailes 1996). Rather, for our pur- 



pose it was sufficient to simulate the mass-exchange episode by 
subjecting a RGB mod el to a sufficientl y larg e mass loss rate 
(Iben & Tutukov |1986| , Castellani et al. 1994 ). We emphasize 
that, once mass loss has been turned off, the model's further 
evolution does not depend on the details of the previous mass 
loss episode. The evolutionary speed across the Hertzsprung- 
Russell diagram is rather controlled by the model's structural 
changes caused by the actual envelope-mass reduction due to 
hydrogen burning and mass loss. 

In order to get realistic starting models, we calculated a 
1 M© sequence from the pre-main sequence phase up to the 
tip of the RGB. Along th e RGB we applied mass-loss rates 
Mr acco rding to Reimers ( 1975 ) with rj = 0.5 as in Maeder & 
Meynet ( 1989 ). At appropriate positions high mass loss rates, 
Mhi g h, were invoked in order to get models of desired final 
mass, M. Since in the actual situations Roche-lobe overflow is 
assumed to occur on a nuclear time scale, the maximum applied 
mass loss rate was chosen in such a way as not to destroy the 
model's thermal equilibrium. Accordingly, MhXgh varied from 
Mhigh ~ 10~ 9 M Q yr _1 for M « 0.15 M to about M h i g h 
w 10~ 6 M Q yr" 1 for M « 0.4 M . Below a critical envelope 
mass, M£"j, sufficient densities and temperatures to continue 
hydrogen burning can only be maintained if the envelope con- 
tracts. Then, the model starts to leave the RGB. At this point 
of evolution mass loss was virtually switched off by decreasing 
M over a short transition period until M = Mr was reached. 
In general we chose T e g = 5000 K to be the point where M 
should have reached the Reimers rate with rj = 0.5. For mod- 
els with M <, 0.2 M Q we had to increase this temperature to 
T c ff = 10000 K because models on the lower part of the RGB 
have temperatures still too close to T e g = 5000 K. 

The critical envelope mass, M™^, which marks the transi- 
tion between expansion and contraction of the star's envelope, 
depends on the mass of the the hydrogen-exhausted core, M c , 
in a way similar to what is found for post-AGB remnants: the 
larger the core, the smaller the residual envelope mass. When 
this value is reached, Roche-lobe contact will be shut off, and 
the remnant will continue its evolution with L as Lfjyd ~ const 
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and R < towards larger T c s . The models' evolutionary speed 
is determined by the evolution of the envelope mass which is 
reduced both by hydrogen burning, i.e. the core growth rate 
M c , and the mass loss rate, i.e. Mr in the presented models. 
Since M c / Mr ~ 1 /R mass loss becomes rapidly unimportant 
at least for T e g > 10000 K. The envelope mass at the turn- 
around point at maximum effective temperature is independent 
of the post-RGB mass loss. 

Our procedure is adequate for obtaining reliable starting 
models for He- WD tracks. Due to their evolutionary history the 
internal structure of these initial models is consistent with the 
one expected from binary evolution (see references in Sect. |]). 
This concerns the envelope masses as well. As will be shown in 
the following sections the latter has consequences for the mass- 
radius relation. Furthermore, the evolutionary envelope masses 
give rise to ongoing hydrogen shell burning during the white 
dwarf evolution. This greatly prolongs the cooling times of the 
He-WDs. 



Table 1. Total remnant mass, M, mass of the hydrogen- 
exhausted core, M c , total mass of the outer hydrogen lay- 
ers ("thickness"), Mh, envelope mass, M cnv , and helium sur- 
face abundance by mass fraction, Y, at T c g = 5000 K for 
M > O.2M and at 10000 K for M < 0.2M Q after the end of 
RGB evolution 



M/Mq 


Af c /M 


M H 
1O- 3 M 


Mcnv 

10- a M Q 


Y 


0.179 


0.1552 


25.550 


48.289 


0.464 


0.195 


0.1782 


16.766 


30.768 


0.462 


0.234 


0.2220 


8.118 


13.098 


0.354 


0.259 


0.2524 


4.771 


7.232 


0.312 


0.300 


0.2960 


3.189 


4.746 


0.301 


0.331 


0.3281 


2.509 


3.744 


0.301 


0.414 


0.4116 


1.446 


2.175 


0.301 



3. Results 

We calculated seven evolutionary sequences starting at differ- 
ent locations along the RGB of the 1 M Q sequence. They are 
listed in Table [IJ, where we give the total amount of hydrogen 
contained in the envelope, the 'thickness' Mh, the total enve- 
lope mass, M cnv , and the helium surface abundance by mass 
fraction, Y. Note, that Mh ~ 0.7 M cnv except for the least 
massive remnants. For M ^ O.3M0 mass loss has uncovered 
layers where hydrogen burning took place during the main se- 
quence phase. Consequently, the helium surface abundance of 
those models amounts up to Y ~ 0.46. M env , and Mh, are 
reduced (and the core mass correspondingly increased) by con- 
tinuing hydrogen burning. Mass loss does scarcely reduce the 
total mass (see Sect. 2.1). 

At T cff = 5000 K and 10000 K, resp., after the end of 
the RGB evolution the envelope masses of the Pop. I He- 
WD models are quite large, varying between 2 • 1O~ 3 M and 



50 • 1O~ 3 M (see Table [l]). It's noteworthy that M cnv would be 
even larger for lower metallicity (Castellani et al. 1994). 

Figure [l] shows the complete evolutionary tracks of all se- 
quences down to about T e g » 4000 K and L w 10 -4 L Q . 
The corresponding data[] are given in Table 5-11. Only for two 
model sequences, 0.234 M and 0.259 M , we found thermal 
instabilities of the hydrogen-burning shell when CNO burning 
ceases (e.g. Kippenhahn et al 



1968 



Iben & Tutukov 



1986) 



The hook-like inversions on the cooling tracks for M = 0.300 
and M = 0.331 M are also due to the onset of unstable burn- 
ing. Full-scale instabilities are, however, avoided because the 
shell regions cool off too effectively. A detailed discussion of 
the properties of these hydrogen flashes is deferred to a sepa- 
rate paper. Here we only note that we adjusted our numerical 
time steps properly in order to resolve these flash phases with 
satisfactory accuracy. It is interesting to remark that no instabil- 
ities occurred below M 0.20 M . The existence of a lower 
limit for the occurrenc e of h ydrogen flashes, which has been 
predicted by Webbink ( 1975 ), has important ramifications for 
any discussions of the evolution of low-mass white dwarf com- 
panions of millisecond pulsars, as will be shown below. 

Before discussing the properties of our helium white dwarf 
models, the influence of different input physics should briefly 
be reported. Concerning treatment of nuclear reactions we 
found no significant differences between tracks (or mass-radius 
relations) using the nuclear network and those only calculated 
with equilibrium rates for the four most important elements (H, 
He, C, O). 

The comparison between tracks calculated with OPAL and 
the older CS opacities yielded, in short, the following differ- 
ences: 

- A well known shift of the RGB locus towards higher effec- 
tive temperatures for CS opacities. 

- Slightly larger turn-around temperatures for the CS post- 
RGB tracks. For example, a CS model sequence with M = 
0.412 M reached logT c ff = 4.915, the corresponding 
OPAL model (M = 0.414 M©) only logoff = 4.911. 

- The opacity effects on the locations of the cooling tracks 
as well as on the mass-radius relation and envelope masses 
were found to be small. Accordingly, the opacities seem 
not to be a critical parameter for the determination of the 
mass-radius relation (for given metallicity). 

- White dwarf cooling ages are, as expected, significantly ef- 
fected by the opacities. Models calculated with CS opac- 
ities have shorter cooling ages than those calculated with 
OPAL. The influence, however, is only noticeable in our 
more massive models where hydrogen burning is less dom- 
inant (see next section). 



3.1. White dwarf cooling properties 

For CO white dwarfs it is well known that the energy contribu- 
tion of nuclear burning drops rapidly below the gravothermal 



Tables 5-11 are only available in electronic form at the CDS via 
ftp 130.79.128.5. 
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Fig. 1. HRD with complete evolutionary tracks of RGB remnants with different masses (from top: 0.414, 0.331, 0.300, 0.259, 
0.234, 0.195, 0.179 M©). The long-dashed curve shows the evolutionary track of the 1 Mq star we used for abstracting the 
remnants by mass loss (see Sect. 2.1). The short-dashed loops mark the redward short excursions due to hydrogen shell flashes 
that occurred for the 0.259 and 0.234 M© sequences when the CNO cycle is shut off. 



energy release (compressional heating plus cooling ) when the 
cooling branch is reached (Iben & Tutukov 1984b; Koester & 



Schonberner 1986; Blocker 1995b). In contrast, full evolution 



ary calculations for helium white dwarfs show that hydrogen 
shell burning still plays a significant role on the cooling branch 
leading to hydrogen shell flashes for certain masses (e.g. Kip- 
penhahn et al. 1968, Iben & Tutukov 1986, Castellani & Castel- 
lani 1993). The present calculations illustrate that hydrogen 
shell burning remains the main energy source down to temper- 
atures well below T c ff = 10 000 K on the cooling branch. For 
demonstration, the temperatures in the burning shells are listed 
in Table |. Also given are the corresponding envelope masses. 
The comparison with those at maximum effective temperatures 
illustrates clearly the significance of hydrogen burning. 



Figure || shows the ratio of the hydrogen shell luminosity, 
^Hyd, to the gravothermal luminosity, L g , for our helium-white 
dwarf sequences as a function of T e g. The shell flashes are 
omitted for clarity. Figure || shows the same ratio as a function 
of the cooling age r. 

The left part of Fig. ^| indicates a rapid drop of L^ y ^/L g 
due to the decline of CNO burning right after the models have 
entered their cooling branches. At lower luminosities and tem- 
peratures, Lnyd/L s exceeds unity again, entirely due to pp 
burning (see below). This phase of dominant hydrogen burn- 
ing lasts for about 20 Gyr for the two least massive models (see 
Fig. ||). Even the 0.414 M Q model burns hydrogen for about 
3 Gyr until gravothermal energy release finally resumes. 
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Table 2. Temperatures at the center and the lower and upper boundary of the hydrogen burning shell for our different white 
dwarf models at T e g = 10 000 K on the cooling branch. Tbot and T top are the temperatures (in K) of those layers where the 
energy generation rate for hydrogen burning, en, has dropped to 1 % of the maximum value. M env knee and M env are the envelope 
masses at maximum effective temperature (T ff kncc) and at T cff = 10 000 K on the cooling branch, resp. Tkncc is the surface 
luminosity corresponding to T e g = Teff^nec- 



M/Mq 


log(T c /K) 


log (T bot /K) 


log(T top /K) 


M cnv 
10- 3 Mg 


M onVjkncc 
lO- a M 


l0gT e ff,kncc/K 


log(L k ncc/i0) 


0.179 


7.092 


7.086 


6.796 


1.460 


5.369 


4.1372 


-0.6202 


0.195 


7.070 


7.068 


6.784 


1.184 


5.214 


4.2100 


-0.2262 


0.234 


7.039 


7.029 


6.766 


0.744 


2.836 


4.4430 


0.8708 


0.259 


7.022 


7.011 


6.748 


0.666 


2.593 


4.5168 


1.1283 


0.300 


7.000 


6.988 


6.730 


0.529 


1.652 


4.6038 


1.4364 


0.331 


6.987 


6.972 


6.721 


0.438 


1.216 


4.7053 


1.7950 


0.414 


6.959 


6.934 


6.704 


0.317 


0.612 


4.9124 


2.4873 




Fig. 2. Ratio of hydrogen to gravothermal luminosity as a func- 
tion of T ff for He-WDs of different masses. 



A more detailed picture of the temporal evolution of the 
different luminosity contributions (Figs. 5] and |]) reveals that at 
higher luminosities almost the entire energy production comes 
from shell CNO burning. At about r w 2 • 10 s yr (0.195 M©) 
and t sa 1 • 10 5 yr (0.414 M ) the cooling branch is reached, 
and L, Tnyd and Tcno start to decrease, slowly in the low- 
mass model, but very rapidly in the 0.414 M model. In the 
less massive models hydrogen burning via the pp chains re- 
mains the predominant luminosity contribution for more than 
10 Gyr until finally gravothermal energy release starts to take 
over (see Fig. Q). Since unstable burning is restricted to a cer- 
tain mass range the feature of continuing pp burning of low 
mass He-WDs (Al ^ O.2M ) is particularly independent of 
the details in modelling the phase of hydrogen shell flashes 
(e.g. the mass loss during the reexpansion phase due to Roche 
lobe overflow). For more massive models, pp burning stays im- 
portant only for a few Gyr (Fig. ||). 

The more gradual transition from CNO burning to pp burn- 
ing is typical for less massive He-WDs. The more massive 
models show a rapid, more step-like shaped transition, very 



Fig. 3. Ratio of hydrogen to gravothermal luminosity as a func- 
tion of the cooling age r. Only the lower parts of the respective 
cooling curves (T eff < 15000 . . . 18000 K) are plotted. The 
ages are counted from T cff = 10 000 K (M < 0.2 M Q ) and 
T cff = 5 000 K (M > 0.2 M ) on the (horizontal) post-RGB 
branch shown in Fig. 1. 



similar to the situation found in CO-WDs (Blocker 1995b). 
Figures ^| and Q together with Fig. ||, also demonstrate that 
neither the gravothermal contribution nor the energy loss via 
neutrinos plays a significant role in the cooling history of He- 
WDs with M <, 0.2 M for ages r <, 10 Gyr. 

Figure ^illustrates how the the ratio Tcno /L pp varies dur- 
ing the course of evolution. The hooks indicate the onset of un- 
stable CNO burning, the flash loops are omitted for clarity. Fig- 
ure H elucidates also that pp burning becomes the dominant en- 
ergy contribution for He-WDs below T cS ~ 35 000 K for M = 
0.414 M Q and below T cff » 18 000 K for M = 0.234 M . 
This coincides with the change from CNO to pp burning as can 
be seen from the minima in Fig. |[ 

Summarizing, nuclear burning remains an important, if not 
dominant, energy source for helium-white dwarfs, even for 
temperatures of T e g = 10 000 K and below. This is in partic- 
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ular true for objects with M < 0.2 M©, like the companion of 
PSR J1012+5307. For example, Fig. | shows that for the WD 
companion of PSR J1012+5307 with M « 0.2 M the ratio 
Lnyd/Lg is at least 50. The assumption £ nU cicar ~ usually 
made in contraction models thus appears not to be appropriate 
forHe-WDs. (e.g. Chin & Stothers[l97l|; Althaus & Benvenuto 



1997) 



H oweve r, Hansen & Phinney ( 1998a ) and Benvenuto & Al- 
thaus ( 1998 ) did consider nuclear burning but found it of little 
importance. They did not model the evolution of the He-WD 
progenitor and had to make assump tions a bout the hydrogen 
envelope mass. Hansen & Phinney (1998a) selected an enve- 
lope mass of Mh = 3 • 10~ 4 M Q taken at some stage dur- 
ing the end of the second hydrogen shell flash from a 0.3 M 
model sequence calculated by Iben & Tutukov (1986). At this 
phase the envelope mass has already been considerably re- 
duced by the hydrogen burning. Additionally, the Iben & Tu- 
tukov (1986) model had suffered from severe mass loss due to 
Roche lobe overflow during the expansion back to the RGB 
domain. Note that in models with M w O.3M and such a 
small hydrogen envelope mass nuclear burning indeed becomes 
less important. However, as mentioned above, shell flashes do 
not occur for M Ss 0.20 M which is the relevant mass range 
for the PSR J1012+5307 companion. The evolutionary mod- 
els show that this object should have a hydrogen envelope of 
Mh ~ 10~ 3 M . Thus, the finding by Hansen & Phinney 



( |1998a| ) that nuclear burning is negligible results from their 
choice of the hydrogen envelope mass which is not appropiate 
according to our models. 

Additionally, the evolutionary history plays an important 
role as well since it determines the thermomechanical structure 
on the prevailing part of the cooling branch which is different 
if the pre-He-WD evolution is not considered. For instance, in 
contrast to full evolutionary calculations (e.g. Iben & Tutukov 
1986, present paper) hydrogen burning can be found to be 
unimportant even for thick-envelope models (Mh ~ 10 -3 M 
for M = 0.2 M©) if the pre-He-WD evolution is not consid- 
ered (Benvenuto & Althaus |1998b. 



3.2. Age and mass of the PSR J '1012+5307 'white dwarf 
companion 

In order to estimate the influence of initial conditions on 
age and mass determinations for low-mass white dwarfs, we 
also calculated contraction models (L nuc = 0) with identical 
masses and chemical profiles as our evolutionary models. For 
that purpose we started with a homogeneous main sequence 
model of a given mass whose chemical profile was adapted 
over a series of about 50 models to the chemical structure of 
the corresponding evolutionary m odel a t the beginning of the 
cooling branch (see Blocker et al. 1997 ), but did not allow for 
hydrogen burning. 
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Fig. 4. Different luminosity contributions Li as a function of 
cooling time r for the 0.195 M -sequence. Ages are counted 
from T off = 10 000 K on the horizontal post-RGB branch. 
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Fig. 5. Same as Fig. |J but for the 0.414 M -sequence. Ages 
are counted from T e g = 5 000 K on the horizontal post-RGB 
branch. 



3.2. 1 . Cooling behaviour of helium- white dwarfs 

Figure ^illustrates in a log g, T e g diagram the differences in the 
cooling properties between contraction models and one of our 
evolutionary models (0.195 M ). At the position of the PSR 
J1012+5307 companion, our evolutionary model indicates an 
age of about 6 ± lGyr. In contrast, our contraction model as 
well as the models of Althaus & Benvenuto (1997) give only 
a cooling age of 0.4 Gyr, about a factor of 15 smaller and in 
serious conflict with the pulsar's spin-down age. Note that the 
contraction models agree in the cooling ag e but differ consider- 
ably in the mas s-radi us-relation (see Sect. 3.2.2 ), since Althaus 
& Benvenuto ( 1997 ) used Mh = while we took the same 
Mh as in the corresponding evolutionary model. 

A more detailed view of the cooling behaviour of our evo- 
lutionary models is given in Fig. || which shows the tracks as 
well as the corresponding isochrones in the \ogg - logT c ff 
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Fig. 7. log g — log T ff diagram with our tracks for M — . 1 9 5 M© with and without consideration of nuclear burning (continuous 
and dotted line, see text). Also shown are the tracks of Althaus & Benvenuto (1997) for M = 0.1, 0.15 and 0.2 M Q (AB, dashed 



lines) as well as time marks for different cooling stages as indicated. Additional circles mark cooling ages of our evolutionary 
model of 1 and 10 Gyr, resp. All ages are counted from a pre-white dwarf stage close to T e g = 10 000 K. The filled triangle 
marks the position of the PSR J1012+5307 white dwarf. 



plane. Note that hydrogen burning leads to different slopes of 
the isochrones for CO- and He-WDs. Usually, more massive 
white dwarfs cool faster then less massive ones. For He-WDs, 
however, the reverse is true. Due to the increasing importance 
of nuclear burning with decreasing mass, the less massive ones 
are considerably older than the more massive ones, at a given 
temperature below about 10000 K. Hence, the isochrones of 
He-WDs run almost perpendicular to those of CO-WDs, and 
an extrapolation of cooling times of CO-WDs into the He-WD 
regime would give significantly lower cooling ages, close to 
those of the contraction models. 

As already mentioned above, our evolutionary models pre- 
dict a cooling age of 6 ± lGyr for the white dwarf compan- 
ion of PSR J1012+5307, which is in good agreement with the 
estimated pulsar's spin-down age of 7 Gyr, and also with age 
estimates of Alberts et al. (1990 who also included nuclear 
burning. 

Thus, evolutionary models for low-mass white dwarfs 
which start from explicit model calculations of the progenitors 
evolution seem to indicate that neither a modification of the 



theory of field decay as proposed by Burderi et al. (1996) nor 
the consideration of initial spin periods close to the present one 
are necessary to resolve the apparent age discrepancy between 
the two components of PSR J1012+5307. 

3.2.2. Mass-radius relations of helium-white dwarfs 

The mass-radius relation according to our models is shown in 
Fig. ^ and also listed in Table ^| For intermediate temperatures 
some distortions can be seen where hydrogen flashes with the 
associated high burning rates lead to a substantial reduction of 
the envelope mass. A more thorough discussion will be pre- 
sented in a forthcoming paper. Important for the interpretation 
of observations is the remarkable temperature dependence of 
the radii even at the low temperature end. 

The structural differences between evolutionary and purely 
contracting white dwarf models influence the evolutionary 
tracks as well. For selected masses the tracks of evolutionary 
and contracting white dwarfs with identical chemical struc- 
tures are plotted in Fig. [k]. The long-dashed line in the fig- 



9 



Table 3. Radii for He-WDs for different masses and effective temperatures (see also Fig. H). The radii are given in solar units. 



M/Mq 





0.179 


0.195 


0.234 


0.259 


0.300 


0.331 


0.414 


Teff/K 


40000 










0.0972 


0.0525 


0.0330 


35 000 










0.0626 


0.0460 


0.0301 


30000 








0.0698 


0.0512 


0.0407 


0.0267 


25 000 






0.0698 


0.0523 


0.0433 


0.0353 


0.0226 


20000 






0.0505 


0.0424 


0.0284 


0.0248 


0.0196 


15 000 




0.0538 


0.0322 


0.0286 


0.0244 


0.0220 


0.0182 


10000 


0.0384 


0.0332 


0.0264 


0.0244 


0.0217 


0.0201 


0.0170 


5 000 


0.0261 


0.0246 


0.0219 


0.0208 


0.0192 


0.0181 


0.0159 




Fig. 6. Contribution of the two hydrogen burning reaction 
chains: ratio of luminosity due to the pp-chain (L pp ) and due to 
the CNO-cycle (Lcno) as a function of T e g on the WD cooling 
branch. The shell flashes are omitted for clarity. 



ure connects roughly the loci where the tracks of evolutionary 
and contraction models do converge. In the regime above this 
line the thermo-mechanical structure of white dwarfs does de- 
pend on their (assumed) evolutionary history. This has impor- 
tant consequences for mass determinations from spectroscopy 
(see Blocker et al. 



1997). If an object's position lies above 



the line its mass will be overestimated by using simple con- 
traction models with the same envelope masses as the corre- 
sponding evolurionary models. This fact is especially impor- 
tant for helium-white dwarfs with their low masses: The tracks 
of evolutionary and contraction models in the log g — log T e g- 
diagram merge only very late at rather low effective tempera- 
tures. For example, the merging of both 0.2 M Q tracks occurs 
atlogT off «3.5! 

Although in the lower right part of Fig. 
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the thermo- 
mechanical structures of evolutionary and contracting mod- 
els agree, the evolutionary speeds may be very different be- 



Fig. 8. log g — log T e g diagram with our evolutionary tracks 
for He-WDs (M = 0.179, 0.195, 0.234, 0.259, 0.300, 0.331, 
0.414M© and for CO-WDs (AT = 0.524, 0.605, 0.696, 0.836, 
0.940M Q ; Blocker 1995b). The masses increase from top to 
bottom. Also shown is the locus of the companion of the mil- 
lisecond pulsar PSR J 101 2+5307 as given by van Kerkwijk 
et al. (1996). Isochrones are shown for 3 • 10 8 ,1 • 10 9 ,3 ■ 
10 9 , and 10 10 yr (from left to right, last isochrone only for He- 
WDs). The isochrones for the He-WDs are turned over and 
shifted to the left because hydrogen burning slows down the 
evolution so much. The dotted isochrone of 6 • 10 9 yr fits the 
position of the PSR J1012+5307 companion. 



cause hydrogen burning slows down the 'cooling' for the lower 
masses. 

The possible differences in structure between evolutionary 
and contracting white dwarf models translate also into differ- 
ences in the mass-radius relations, as demonstrated in Fig. pT[ 
One can again see that according to contraction models the ob- 
ject's mass for given gravity is overestimated. However, the 
larger hydrogen envelope masses which result in our models 
from their evolutionary history have the effect to increase the 
masses derived from atmospheric parameters and mass-radius 
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4.5 4 

LogT eff [K] 



Fig. 9. Mass-radius relation for He-WDs at different effective 
temperatures, calculated from our evolutionary models from 
T ofr = 40000 K (upper line) down to T off = 5000 K (lower 
line) in steps of AT c ff = 5000 K (Table ||). For comparison also 
the T — K relation of Hamada & Salpeter (1961) for pure 
helium models (dashed line) is given. Along the dotted lines 
gravity is constant, with (from above) \ogg = 6.0, 6.5, 7.0 and 
7.5). Temperatures larger than 10000 K are not reached by 
RGB remnants of too low a mass (see Fig. [j]). The irregular- 
ities around T e s = 20000 K are due to thermal instabilities in 
the hydrogen burning shell as explained in the text. 



relation. Thus the two effects of contraction vs. evolutionary 
models (with regard to the treatment of L nuc = 0) and smaller, 
mass independent envelope masses vs. evolutionary (thicker, 
mass dependent) envelope masses compensate each other to 
a certain extent. This can be seen from Fig. |l2| which gives 
mass-radius relations of different origin in order to compare 
the mass determination of the PSR J1012+5307 companion. 
The deviations from the T — K relation of Hamada & 



Salpeter ( 1961 ) for a pure helium composition demonstrate 
the well known importance of finite temperature effects in the 
equation of state, especially for the low-mass WDs. For a sur- 
face gravity log g = 6.75 of the pulsar's companion that re- 
lation yields M ps 0.12 M Q which is close to the lower limit 
M • sin i — 0.11 Mq. The differences between our mass-radius 



relation and those of van Kerkwijk et al. (1996) and Althaus & 
Benvenuto ( 1 997 ) are mainly due to the larger envelope masses 
(i.e. thicker hydrogen layers) predicted from our evolutionary 
models. For example, van Kerkwijk et al. (1996) interpolated 
between the Hamada & Salpeter relation based on pure helium 
cores (Mh = 0) and that based on models from Wood (1994) 
with 'thick', but constant hydrogen layers (Mh = 10 -4 Af»). 
Our evolutionary calculations, however, revealed a large de- 
pendence of envelope thicknesses on total masses: Mh drops 
from ps 10~ 3 M© for M = 0.2 M to M H « 10~ 4 M for 
a C/O WD with M = 0.6 Mq, the heaviest He-WDs with 
M « 0.45 M Q havi ng ab out M H = lO" 375 M Q (see Ta- 
ble 3. Blocker et al. |l997| ). Note that in contrast to Fi gs. |io| 
and [ll] the contraction models of Althaus & Benvenuto (1997) 



Fig. 10. log g — log r c ff-diagram with WD tracks according to 
evolutionary models (i.e. with consideration of nuclear burn- 
ing, solid curves) and contraction models (i.e. no nuclear burn- 
ing, dotted curves) with the same masses and chemical struc- 
tures for the masses as indicated by the labels. The evolution- 
ary tracks are either from this study or from Blocker ( 1995b| ) 
for the larger masses. The contraction models are from Blocker 
et al. (1997). The dashed line divides the diagram roughly in 



two parts: in the lower right part contraction and evolutionary 
models have virtually the same thermo-mechanical structure. 
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Fig. 11. Mass-radius relation derived from evolutionary 
(crosses) and contraction models (squares) with identical 
chemical profiles at T cfi = 10 000 K and 20000 K. Also given 
are the T = K relation of Hamada & Salpeter ( |l96l[ ) (dots) 
and lines of constant surface gravities (dashed, from above: 
log g = 6.5,7.0,7.5). 



give smaller white dwarf masses compared to our evolutionary 
models due to the much smaller envelope masses (in this case 
M H = 0). 

As can be seen from Fig. |T^, our mass-radius relation 
yields a mass of M = 0.19 ± 0.02 M Q for the white dwarf 
companion of PSR J1012+5307, using \ogg = 6.75. This is 
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Fig. 12. Mass-radius relation at T e g — 8550 K, relevant for 
the mass determination of the PSR J1012+5307) companion, 
from different authors: AB = Althaus & Benvenuto (1997), K 



= van Kerkwijk et al. ( |1996| ); HS = Hamada & Salpeter (|1961j). 
Lines for constant surface gravity are also given, starting from 
logg = 6.75 (upper line, relevant for PSR J1012+5307) to 
log<? = 7.75 in steps of A log g — 0.25. 



larger than the value quoted by van Kerkwijk et al. ( 199f ) and 



Callanan et al. (1998) but within the mass range A/h -wd = 
0.13 .. . 0.21 M given by Hansen & Phinney ( |1998b| >. Simi- 
larly, A f — 0.16 M follows from the Althaus & Benvenuto 



( |1997[ ) relati on. With the spectroscopic analysis from Callanan 
et al. ( |1998| ) a mass of M = 0.15 ± 0.02 M follows by ex- 
trapolating our mass-radius relation. 

With the mass ratio M pu i sar /A/H C -WD ~ 9.5 ± 0.3 (van 
Kerkwijk 1998, priv. comm.) the pulsar's mass can be es- 
timated: A/hc-wd = 0.15 M gives Af pu ] sar « 1.43 ± 
0.25 M , close to the canonical pulsar mass of w 1.5 M Q . 
On the other hand, A/hc-wd = 0.19 M Q gives Af pu i sar w 
1.81 ± 0.25 M Q . Slightly increased values for M pu i sar are 
obtained if the mass ratio M pu i sar /MH e -WD ~ 10.5 ± 0.5 



from Callanan et al. ( |1998[ ) is applied. A summary of the re- 
sults for the mass determination of the components of the PSR 
J1012+5307 system is given in Table ^. 

4. Summary and conclusions 

We have computed a representative grid of evolutionary white 
dwarf models of low mass for solar composition whose struc- 
ture is consistent with their earlier history of being remnants of 
the first giant branch. Though such models have already been 
presented in the literature for isolated cases, our calculations 
are the first that cover a whole range of masses. These models 
allowed then to draw several important conclusions: 

- We found thermal instabilities due to hydrogen burning 
via the CNO cycle in a geometrically thin shell (hydro- 
gen shell flashes) at the beginning of the cooling branch of 
our evolutionary tracks. It turned out, however, that these 



Table 4. Spectroscopic results, i.e. log<? and T e g, and de- 
rived masses, Af Hc -WD, for the PSR J1012+5307 He-WD. 
The lower part gives the He-WD mass from the present work 
(A/hc-wd,pw) derived from the respective spectroscopic anal- 
ysis. Afp U i sar . k is the pulsar mass based on the mass ra- 
tio Afpuisar/A/Hc-WD ~ 9.5 ± 0.3 (van Kerkwijk et al. 
1998, priv. comm.) and Af Ho -WD,PW- Afp U isar,c is th e pul - 
sar mass based on the mass ratio from Callanan et al. ( 1998[ ) 
(Afpuisar /A/hc-wd ~ 10.5 ± 0.5) and Af Ho -WD,PW- 



Callanan et al. 



(1998) 



van Kerkwijk et 
al. (1996) 



logg[cm/s 2 ] 
T efl [K] 
Mh6-wd[Mq] 



6.34 ±0.20 
8670 ± 300 
0.16 ±0.02 



6.75 ±0.07 
8550 ± 25 
0.16 ±0.02 



Mhe-wd,pw[Mq] 

Afpuisar, K [M ] 

M pu i sar , c [M Q ] 



0.15 ±0.02 
1.43 ±0.25 
1.59 ±0.30 



0.19 ±0.02 
1.81 ±0.25 
2.00 ± 0.30 



shell flashes do occur only within a limited mass range, 
Af = 0.21 . . .0.30 M , independently from the detailed 
input physics used. 
- The cooling speed is reduced by ongoing hydrogen burn- 
ing at the bottom of the envelope via the pp chains. This 
effect is the more pronounced the less massive the models 
are and may slow down the cooling down to 10 000 K by up 
to a factor of about 40. This is important for age determina- 
tions of binary systems containing helium-white dwarfs by 
using the cooling age of the white dwarf components. The 
employment of inconsistent white dwarf models may give 
grossly wrong results in the sense that ages are underesti- 
mated. 

We re-determined the age of the white dwarf component in 
the PSR J1012+5307 system with our evolutionary models 
and found an age of 6 ± lGyr, 15 times larger than that 
estimated from contraction models, but in excellent agree- 
ment with the pulsar's spin-down age of about 7 Gyr. We 
confirm the results of Alberts et al. (1996 1 who also found 



7 Gyr based on their white dwarf models. 
- Our evolutionary models predict an inverse correlation of 
envelope mass with total mass, which is a continuation of 
the correlation already apparent for the more massive white 



dwarfs with carbon-oxygen cores (Blocker et al. 1997). At 
the low mass end, the evolutionary envelope masses are at 
least ten times larger than the often ad hoc assumed values 
for contraction models. Translated into a mass-radius rela- 
tion, these larger envelope masses lead to correspondingly 
higher mass estimates for any given g, T c g, and especially 
for the white dwarf companion of PSR J1012+5307 a mass 
of Af = 0.19 M follows, instead of Af = 0.16 M Q by 
using existing contraction models. 
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